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Abstract 

In this paper, we establish global existence of smooth solutions for the Cauchy prob- 
lem of the critical nonlinear wave equation with time dependent variable coefficients in 
three space dimensions 

du(t> ~ d Xi {g l3 {t, x)d X] 4>) + 4? = 0, K ( xl^ 

where [gij(t,x)j is a regular function valued in the spacetime of 3 x 3 positive definite 
matrix and (jf^it, #)) its inverse matrix. Here and in the sequence, a repeated sum on 
an index in lower and upper position is never indicated. In the constant coefficients case, 
the result of global existence is due to Grillakis [1] ; and in the time- independent variable 
coefficients case, the result of global existence and regularity is due to Ibrahim and 
Majdoub [6]. The key point of our proofs is to show that the energy cannot concentrate 
at any point. For that purpose, following Christodoulou and Klainerman [3], we use a 
null frame associated to an optical function to construct a geometric multiplier similar to 
the well-known Morawetz multiplier. Then we use comparison theorem originated from 
Riemannian Geometry to estimate the error terms. Finally, using Strichartz inequality 
due to [T2] as Ibrahim and Majdoub [6], we obtain global existence. 

Keywords: time dependent variable coefficients, critical nonlinearity, null frame, 
comparison theorem. 

1 Introduction 

In this work we consider global existence of smooth solutions of the Cauchy problem 

d a <p-d Xi (g i ^t,x)d Xj (f) +0 5 = 0, R t xR 3 x , 
^ ,rc) = / 1 (x)GC oo (R3), Uto,x) = f 2 (x)€C^(Rl), 

here {g tJ (i, x)} 3 =1 denotes a matrix valued smooth function of the variables (t,x) G R x 
R 3 , which takes values in the real, symmetric, 3x3 matrices, such that for some C > 0, 

C\ti\ 2 < !r'U..r)ZiZj < C- 1 ^ 2 , V^el 3 , (t,x) elxl 3 . (1.2) 
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Obviously it is a critical wave equation on a curved spacetime. First let us survey existence 
and regularity results for critical nonlinear wave equations briefly. If g 13 = S tJ , which denotes 
the Kronecker delta function, we say problem is of constant coefficients. In the case of 
critical nonlinear wave equation with constant coefficients, a wealth of results are available in 
the literature. For cauchy problem, global existence of C 2 -solutions in dimension n = 3 was 
first obtained by Rauch [9], assuming the initial energy to be small. In 1988, also for "large" 
data global C 2 -solutions in dimension n = 3 were shown to exist by Struwe [T4] in the radially 
symmetric case. Grillakis [4] in 1990 was able to remove the latter symmetry assumption 
and obtained the same result. Not much later, Kapitanskii [7] estiblished the existence 
of a unique, partially regular solution for all dimensions. Combining Strichartz inequality 
and Morawetz estimates, Grillakis [5] in 1992 established global existence and regularity 
for dimensions 3 < n < 5 and announced the corresponding results in the radial caes for 
dimensions n < 7. Then Shatah and Struwe [10] obtained global existence and regularity 
for dimensions 3 < n < 7. They also proved the global well-posedness in the energy space 
in [TT] 1994. For the critical exterior problem in dimension 3, Smith and Sogge [13] in 1995 
proved global existence of smooth solutions. In 2008, Burq et all [2 J obtained the same result 
in 3-D bounded domain. 

For the critical Cauchy problem with time-independent variable coefficients, Ibrahim 
and Majdoub [6] in 2003 studied the existence of both global smooth for dimensions 3 < 
n < 6 and Shatah-Struwe's solutions for dimensions n > 3. Recently, we have showed 
global existence and regularity in [T7] for the critical exterior problem with time-independent 
variable coefficients in dimension n = 3. 

In this paper we are interested in the critical case with coefficients depending on the 
time and space variables. Our result concerns global existence and regularity, showed as 
follow: 

Theorem 1.1. Problem (jl.ip admits a unique global solution <fi E C°°(M x M 3 ). 

The demonstration of theorem 1.1 is done by contradiction, showing is uniformly 
bounded. For that purpose, the key step is to show the non-concentration of the L 6 part of 
the energy (and hence the energy), and to do this the idea is to work in geodesic cone just 
like light cone in constant coefficients case. Thus we have 

Lemma 1.2. (Non-concentration lemma) If G ([£{)) 0) x M 3 ) solves (|1.1|) . then 



where Q(t) is the intersection of t time slice with backward solid characteristic cone from 
origin. 

In order to prove the non-concentration lemma, in the constant coefficients caes the 
Morawetz multiplier tdt + rd r + 1 is used, where r = \x\\ while in the time-independent 
variable coefficients case the geometric multiplier tdt + pVp + 1 is used instead, where p is 
the associated distance function. The time-dependent variable coefficients case considered in 
this work is much more difficult, and the simple minded generalization to use multiplier tdt + 
(u — t)V(u — t) + 1 will not work, where u is an optical function (close to t + \x\ )■ Following 
Christodoulou and Klainerman [3] we use a null frame. However, the emphasis in their work 
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is the asymptotic behavior of the null frame at infinity, and here we emphasize its asymptotic 
behavior locally at a possible blow up point. We derive the asymptotic behavior of the null 
frame by using comparison theorem originated from Riemainnian geometry. 

To prove our result, we also need Strichartz inequality, stated as 
Lemma 1.3. (Strichartz inequality) Assuming g l3 {t,x) satisfies the conditions of the 
introduction, (f) solves the Cauchy problem as follow in the half open strip [t ,0) x M 3 : 

{ d tt <j>-d Xi (g ij (t,x)d Xj cj)) =F(t,x), 

\ <Kto,x) = h(x) e C?(R 3 X ), Mto,x) = f 2 (x) e C °°(R|), 



then we have 



2q 



Lf*Ll([to, 0)x 



< Cdl/lllff^R 3 ) + ||/2||l 2 (R3) + \\ F \\L\Ll([t , 0)x R3)) 



Q<q<oo. (1.3) 



For the proof see Smith [12J. 

Then combining these two lemmas we can establish the uniform bounds on the local 
solution (f), which implies our result, this step is completely parallel to Ibrahim and Majdoub 
[6] and we omit it. 

Our results can extend to a more general variable coefficient second order partial differ- 
ential equation with operator as follow: 

L = d 2 t + 2b% x)d 2 ti - a ij (t, x)df j + Li, a ij = a ji , 

where all coefficients are real and C°°, and L\ is a first order operator. In general, we can 
get rid of cross terms (that is, terms like 6*<9^ ) by the following procedure: let us write 
(with new first order terms ) 

L = (8 t + b%) 2 - (a^ + 6^)flg + L' x . 
If, in the region under consideration, we can perform a change of variables 

Xi = (pi(t,x), ■■■ ,X n = tp n (t,x),T = t, 

and set 

at oxi 

in such a way that the vector field dt + b % di becomes 8t, then the operator L takes the form 

L = %-&8$ CiXj +L 1 , 

for some new coefficients a 13 and lower order terms L\. 

As an application of our result, we can consider the critical wave equation in the 
Schwarzschild spacetime (A4,g) with parameter M > 0, where g is the Schwarzschild metric 
whose line element is 



ds 



2 



2M\ , n / 2M\-1 , 9 , n 

1 )dt 2 +(l J dr 2 + r 2 du 2 , 

r J \ r ) 
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where dui 2 is the measure on the sphere s 2 . While the singularity at r = is a true metric 
singularity, we note that the apparent singularity at r = 2M is merely a coordinate singu- 
larity. Indeed, define the Regge- Wheeler tortoise coordinate r* by 

r* = r + 2M log(r - 2M) - 3M - 2M log M, 

and set v = t + r*. Then in the (r*,t,u)) coordinates the Schwarzschild metric g is expressed 
in the form 

ds 2 = - (l - - — ) dv 2 + 2dvdr + r 2 dw 2 . 

Let £ be an arbitrary Cauchy surface for the (maximally extended) Schwarzschild space- 
time (Ai,g) stated as above and consider the Cauchy problem of the wave equation 

|^-* 5 = o, (l4) 

for this problem Marzuola et all [8] proved global existence and uniqueness of finite energy 
solution under the assumption of small initial energy, and according to our result we can 
remove the small energy assumption, that is 

Theorem 1.4. For smooth initial data prescribed on E, equation (|1.4|) admits a unique 
global smooth solution in the (r*,i,w) coordinates. 

Now we sketch the plan of this article. In the next section we recall some geometric 
concepts which are necessary for our proofs. Section 3 is devoted to the proof of lemma 1.2: 
the fundamental lemma expressing the non-concentration of L 6 part of the energy. 

Finally, we remark that although our non-concentration lemma is stated only in dimen- 
sion n = 3, the proof works in any dimension for the critical wave equations. 

In this paper, the letter C denotes a constant which may change from one to the other. 



2 Null frame 

Let {gij(t, x)}f - =1 denotes the inverse matrix of {g^it, £c)}? J=1 , and consider the spilt 
metric g = —dt 2 + gij(t,x)dx t dx : ' = g a pdx a dx^ on R^ t (close to the Minkowski metric). So 
we will work in the spacetime, a 4-dimensional manifold M. Local coordinates on M are 
denoted by x a ,a = 1,2,3,4. The convention is used that Latin indices run from 1 to 3 
while Greek indices relate to the spacetime manifold M and run from 1 to 4. The index 
4 corresponds to the time coordinate, while (xi,X2,xs) are the spatial coordinates. The 
corresponding partial derivatives are d a = We introduce an optical function u (close 
to t + |x| )for g: a C 1 function which satisfies the eikonal equation 

g aP d a udpu = g a pd a ud p u = (Vn, Vn) = |Vn| 2 = 0, 

(2.1) 

u(t,0) = t, 
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where ( , ) denotes the inner product about the given metric. In PDE terms, this means 
that the level surfaces {u = C} are characteristic surfaces for any operator with principal 
symbol <7 Q/3 Ca?/3- From this construction it is easy to see that the first order derivatives 
of u are locally bounded. 
Then we set 

L = -Vu = (d t u)d t - {g ij d l u)d ] = m~\d t + N), 

T dt , (sm^dj ( 2 - 2 ) 

L = Wu + ^MT = m{dt ~ Nl 

where V is the gradient about the given metric, m = (dtu)^ 1 , N = — — -§0^- = — (mg^diUjdj. 
It is easily to see that they are close to dt — d r and dt + d r respectively. And Dl_L = 0, show- 
ing that a integral curve of L is a geodesic. This follows from the symmetry of the Hessian, 
since for any vector field X, we have 

< D L L,X >=- < DlVu,X >= - < D x Vu,L > 

=< D X L,L >= -X <L,L >= -X < Vu,Vu >= 0. 

So the integral curves of the field L generate a backward geodesic cone with vertices on the t- 
axis. Using the coordinate t, we define the foliation J2t ± = {( x ^)^ = an d using u, we 
define the foliation by nonstandard 2-spheres as 

= {{x,t),t = h,u(x,t) = nj. 

Since Vu is orthogonal to {u = u x } and dt is orthogonal to Ylti > ^ e — ^ s a nun vec tor 
orthogonal to the geodesic cone and is an horizontal field orthogonal to <St 1)Ml . Moreover, 

(N,N) = -^—g^dku^dm) = -l-g^dkudm = 1. 

(d t u) {dtu) 2 
Then, if (ei, e<i) form an orthonormal basis on the nonstandard spheres, the frame 

ei,e 2 , e 3 = L = m(d t - N), e 4 = L = m^ 1 (d t + N) 

is a null frame with 

(ei, ei) = (e 2 , e 2 ) = 1, (ei, e 2 ) = 0, 
(ei,L) = (ei,L) = (e 2 ,L) = (e 2 ,L) = 0, 
l(^^) = a,^) = 0, (L,L) = -2. 
We will work in the null frame as above and it requires that we know the vector D a ep, that 
is: the frame coefficients < D a ep,e 1 >. 
We define the frame coefficients by 

< D a L, e b >= x ab = X ba , < D a L, e h >= Xab = Xba, 

< D L L, e a >= 0, < D L L, e a >= 2£ a , 

(2.3) 

< D L L, e a >= 27] a , < D L L, e a >= 2r? a , 
<D L L,L>=0, < D L L,L>= 4w = - < D L L,L >, 
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where a, b = 1, 2. 

If we call k the second fundamental form of Ej by 

fcpT, y) = - < £>x^, Y >, % = 

then k is nothing but the first order derivatives of g and so bounded. By some simple 
computation, we also have 

2f? a = -2fc Ar a , 
2r/ a = -2me a (u t ) + 2% a , 
2£ a = -2m 2 r]^ + 2m 2 k Na , 



2.4) 



Xa fe = -m 2 y,-2m^, 



co = -d t m = . 

(dtu) 2 

For the details, one can read Alinhac jTj. And we are interested in the asymptotic behavior 
of these frame coefficients near the origin, thus we have 

Theorem 2.1. Assuming £ a , r] a , r] , to, x_ ab , Xab are frame coefficients as above, then 



M < c, 

ct 

— < to = 
2 



(hill 

(dtu) 



., - m 2 d u u < ~, 



1 1 

+ ct < x < 



ct, 



t — u - aa t — u 

\v I < -ct, 

'-La' — 

\U<c-ct, 

1 1 — ml I v I < C, 

4 + Ct < XaaU + XbbU + x aa u + X bb u <A~Ct, 



(2 + Cr:)|V0| 2 < £ (XabU + X^^^e^) < (2-Ct)|V0| 2 



(2.5) 
(2.6) 

(2.7) 

(2.8) 
(2.9) 
(2.10) 
(2.11) 

(2-12) 



a,b=l 



where a = 1, 2; c, C are positive constants; t < as we work in the backward geodesic cone 
starting from the origin and u = 2t — u. 

The first inequality is obviously from (|2.4p . and to prove the other inequality of this 
theorem, we need several lemmas stated below. First let us introduce the comparison theo- 
rem. 

Assuming C, D take values in the real, symmetric, (n — 1) x (n — 1) matrices. If for 
any (at,-- - , a n -l), (P\,'" , Ai-l) £ ^ n ~ 1 and «? = I] ft 2 ) we nave 



ai, 



- «n-l 



c 



>(/?!,••• ,/3n-l ) 



V A»-i / 



then we say C >~ D. 

Lemma 2.2. Let gl(n — 1, M) be set of n — 1 order real symmetric matrix, X, : [0, b) 
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gl(n — 1,R). Suppose A : [0,6) — > gl(n — 1,R) satisfies the ordinary differential equation 

d 2 A A 

+ AK = 0, 



ds 2 



dA 



A(0) = 0, —7- (0) = /(i/ie unit matrix), 
ds 



and A : [0, 6) — >• gZ(n — 1,M) satisfies 

d 2 A 



ds 2 



+ AK = 0, 



A(0) = 0, (0) = I (the unit matrix), 
ds 

where s £ [0, 6). Also we assume A, A are invertible in [0, 6) and K -< K, then 

ds ds 



(2.13) 



For the proof see |16| . 

If the metric g = —dt 2 + gij(x)dx l dx^ , depend only on the spatial coordinates, 

then u = t + p is an optical function for g satisfies u(t, 0) = t, where p is the Riemannian 
distance function on the Riemannian manifold (M 3 , gij(x)). The corresponding null frame 
related to u_ is 

ei, e 2 , e 3 = d t + = d t + ^Pidj, e 4 = d t - dp = d t - ]? 3 Pidj, 

where {g tJ ' (x)}fj =1 denotes the inverse matrix of {gij(x)}f =1 . And then u can be compared 
with u through lemma 2.2. 

Let 7 : [0, 6) — > (R^ t , <?) is the integral curve of L = — Vit and we call it null geodesic, 
as < L, L >= 0, then 7 = L = — Vu. Let {ei,e2,e3 = L, e 4 = L} be parallel null frame 
along 7, Ji(s) be the Jacobi field along 7, satisfies Jj(0) = 0, Jj(0) = ej(0),(i = 1,2,3). So 
we have 





where A(s) denotes an invertible matrix valued function of the parameter s £ [0,6). Then 
the Jacobi equation becomes 



cPA 
ds 2 



+ AK = 0, 



where K = (Kij)fj =1 ,Kij =< i?(7, ej)7, > denotes a symmetric matrix of 3 x 3. We then 
easily get 



Hu(ei,ej) = D 2 u(ei,ej) = -(A 1 



dA s 
~ds' 



where denotes the Hessian form of u. And then (|2.3[) yields 



2 Hi \ 



2^11 K12 
K12 K22 ^2 



(2.14) 



(2.15) 
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Correspondingly for optical function u = t + p, we have A, K and 

~ dA 

H^ej) = LPufaej) = D 2 (t + p) = -(A- 1 —)^ = (2.16) 

Note that one assumption of lemma 2.2 is K -< K, but -K33 =< R(L, L)L, L >^ 
0, while ^33 =< R(dt — dp, dt + dp)dt — dp, dt + dp >= 0, so we have to introduce a 
conformally related metric tensor to g to ensure the condition K -< K. Let (gij(x),M. 3 ) be a 
space form with positive constant sectional curvature c. We set then the conformally related 
metric g c = e ct2 g. 

Lemma 2.3. Let (M, g) be a semi-Riemannian manifold of dimension n and let g c = ipg be 
a conformally related metric tensor to g, where <p : M — > (0, 00) is a map. Then 

(1) V = -V, where V and V are the gradients on (M,g) and (M,g c ), respectively. 

(2) For X, Y £ TTM, 



VxY = V X Y + ^-X(p)Y + ^-Y(<p)X - ±g(X, Y)V<p, 
lip lip lip 

c 

where V and V are the Levi-Civita connections of (M,g) and (M,g c ), respectively. 
(3) If / : M -»• R is a map then, for X, Y G TTM, 

H c f (X,Y) =H f (X,Y) - ^-[g(Vp,X)g(Vf,Y) 

+ g(Vf, X)g(Vip, Y) - g(Vp, Vf)g(X, Y)\, 

where Hf and Hj are the Hessian forms of / on (M,g) and (M,g c ), respectively. 
For the proof see [15]. 
From lemma 2.3, we have 

< V(* + P), V(t + P) >g = e ct2 < -^V(i + p), ^V(t + p) >g= 0, 

then i + p is also an optical function for g c . As above, we define A c , K c , Xtij associated to g c . 
It is easily known that for a manifold (M, g^i) with constant curvature c 

R(X, Y, Z, W) 

(2.17) 

= c[g M {X, Z)g M (Y, W) - g M (X, W)g M (Y, Z)] ,VX,Y,W,Ze TTM. 
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Then for the space form (gij(x) 
tation according to (|2.17p gives 



with positive constant sectional curvature c a compu- 





f c 






K = 


c 





■ 




^00 


o) 




*41 ~~ " 


-Hu(ei, 


ei) 


= - < -D ?1 V(i + j5), ei >g= - < £>e x Vt, ei > 






, ei 


) = — v / ccot( v / cp) < ei, ei >^= — y/ccot(y/cp), 


^22 ~~ ~ 


-Hu(e 2 , 


e 2 ) 


= —y/cCOt{y/cp) < e 2 , e 2 >g= — y/c COt ( \/C/5) , 


2£l2 ~~ ~ 


-Hu(ei, 


e 2 ) 


= — v^ccotfVc?) < ei, e 2 >gr= 0, 


^33 ~~ ~ 


-Hu(e3, 


e 3 ) 


= - < Ar 3 V(t + p), e 3 >g=< D dt+d pd t - dp, 


Kl3 ~~ ~ 


-Hu(ei, 


e 3 ) 


= - < D ?3 V(t + p), ex >g=< D dt+d pd t - dp, 


^23 ~~ ~ 


-Hu{e 2 , 


e 3 ) 


= - < De 3 V(t + p),e 2 >g=< D dt+d pd t - dp, 



where V 9 is the gradient on the space form 
conformal change of metric they become 

*2 



. Thanks to lemma 2.3, after the 



2ce 



c 2 t 2 e ct2 







X 



ell 



X 



•r22 



X 



.(■12 



X 



.(■:« 



^cl3 
^c23 






-fll(e 2j e 2 ) 
e 2 ) 

-^|(e 3 , e 3 ) 
-fil(ei, e 3 ) 
-H^2, e 3 ) 



2ce 



cf 2 



c 2 i 2 e rf2 








4ce rf2 
(i^„(ei, ei) - ci) = -Vccot(\/cp) + ct, 



-(Hu(e 2 , e 2 ) - ct) = 
-Hu(ei, e 2 ) = 0, 
-(#g(e 3 , e 3 )-2cf) 
-Hu(ei, e 3 ) = 0, 
-Hu{e 2 , e 3 ) = 0. 



-\/ccot(\/cp) + ct, 
2ct, 



(2.18) 



So we can ensure -fC -< ET C when i 

Similarly, if we let (jjij(x 
ture — c and set the conformally related metric <7 C = e~ ct2 g, then we have 



-7- and c big enough. 

be a space form with negative constant sectional curva- 





f -2ce 


-ct 2 _ c 2 t 2 e ~ct 2 o 





K c = 




-2ce- c * 2 - 


c 2 t 2 e- rf2 




\ 





-4ce" 


2£cii ~~ ~ 




fi) = + = 


= — \/ccoth(\/cjo) 


^c22 ~~ ~ 


-i?f(t 2 , 


f 2 ) = -(#~(t 2 , + = 


= — v / ccoth(v / cpi) 


Kel2 ~ ' 




e 2 ) = --Hg(ei, e 2 ) = 0, 




^c33 ~~ ~ 




f 3 ) = -(fl=(? 3j ? 3 ) + 2ct) 


= -2ct, 


^cl3 ~~ ~ 




e 3 ) = --ffc(ei, e 3 ) = 0, 




2^c23 ~~ ~ 




e 3 ) = -i%(e 2 , e 3 ) = 0. 





(2.19) 
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Again we can ensure K c -< K. 
Using lemma 2.2, we get 

together with (|2.14p . we have 
So 



ds c ds 



(Hu{ei,ej) X (H^ej)). 



4w = - < D L L, L >= Hu(e 3 , e 3 ) < iJ~(e 3 , e 3 ), 
combining (|2.18p and (|2.22[) . we obtain 

ct 



For the same reason we have 



a " 2"' 



(2.20) 

(2.21) 
(2.22) 

(2.23) 
(2.24) 



and the inequality (2.6) of theorem 2.1 follows. 

Combining lemma 2.2, (|2.18p and (|2.19[) . we obtain 

-Vccot(Vcp) + ct = x caa < X aa < X caa = -Vccoth(y/cp) -ct, a = 1, 2, (2.25) 

as we adopt comparison theorem along integral curves of L = —X7u , L = -Vw and L = 
— Vu we set u = u = t + p = u = t + p, so when t, p, p are small (close to 0), we have 



1 1 

+ ct < x < 



t — u 



t — u 



ct, a = 1,2, 



(2.26) 



which is the desired inequality (2.7). 

Using lemma 2.2, together with (|2.14|) . we have 

(^.cij)«j=l ^2Lij) i,j=l (:Kcjj) i,j=l ' 



SO 



° 0^; 



< I 1 o 1| \x l3 




thus 



^(10 1 J &y) 



^cll + 2 2(cl3 + * c3 3 ^ *ll + 2 ^13 + *33 ^ &11 + ^Xcl3 + Xc33- 



(2.27) 
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As from lemma 2.2 we have 



X cii < X# < Xcii, 2 = 1,2,3, 

then ()2.27p means 

Lll + 2 Kcl3 + ^33 ~ Xcll - Xc33 < 2 X 13 < Xdl + 2 X C 13 + X c3 3 ~ X cll ~ X c33 - ( 2 - 2 8) 
Combining (ETT^ . (pETS]) . (gjg)) and (123B1) . a calculation gives 

y<-ct. 

For |jy J we have the same result, and then we obtain the inequality (2.8) in theorem 2.1. 

After that, inequality (2.9) can be easily obtained from (|2.4p . 

To prove the inequality (2.10) of theorem 2.1 we need the following two steps: 


First, show is bounded, where V is the gradient on Euclidean space. From equation 

(|2.ip . we have 

dt(g tJ diudju) = dtg %3 diudjU + 2g lJ 'diudtdjU 
= dtg lJ diudjU + 2g^ diudjU t 
= d t (d t u) 2 
= 2d t ud u u- 

As the first order derivatives of u are bounded, together with (2.6), we have 

\N(u t )\<C. 

Also from (|2.4|) . we get 

me a (u t ) = -rj^ + k Na , 
thanks to inequality (2.8), it implies 

\e a (u t )\ <C, a = 1,2, 

so we finish the first step. 

Second, as the result of the first step 

\l-m\ = I d*~ - 1 1 = \m\\u t (t,x) -u*(*,0)| < Csup| VMtll^l < C\x\. (2.29) 

OtU x 

Set v = t + 5\x\, 5 > 0, then if we choose 5 small enough 

iiv) = L{v) = d t ud t v - g ij diudiv = d t u - 5g ij ' diUpr > 1 - CS > 0, 

\x\ 

while 

7(«) = Lhk) = 0. 

As 7 is a backwards integral curve of L, along the curve 7 we conclude 

u > v = t + 5\x\, 
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thus 

u-t>$|a?|. (2.30) 
Combining ([239]) . ([230]) and (2.7), we have 

11 - m\\x I < Ci, 

which means inequality (2.10). 

Lemma 2.4. Inside the geodesic cone where u < , we have 

\u\<C\t\, \u\<C\t\. (2.31) 
Proof. By ()2.30p . along the integral curve of L starting from the origin, we have 

t < t + 5\x\ < u < 0, (2.32) 

then 

2t < u = 2t - u < t, 

so we complete the proof of lemma 2.4. 
By (pO| . we have 

XaaU + XbbU + X aa U + X^U 



= %aa + ~ 3i) + (! - m2 )(x aa + X^ - 2m(A; aa + k bb )u, 

then (2.7), (2.10) and lemma 2.4 yield the inequality (2.11). 

Now we prove the last inequality of theorem 2.1. Using lemma 2.2 again, we have 

( ei (0) e 2 (0) ) (x^.) e 2 (0 

V o 

/ ei(</>) 
e 2 (</>) 

V o 

/ ei(^) 



^ ( e x (4>) e 2 {4>) J (x^.) 



< ( ei(0) e 2 (0) ) (x c ^.) 



together with <|2.18j> and (j2.19j) . we arrive at 

( - V~ccot(V~cp) + ct) ((e^)) 2 + (ei(<A)) 2 ) 

2 

- Yl X ab e a (<A)efe(</>) 
o,6=l 

< ( - v^coth(V^) - ct) (( ei (0)) 2 + ( ei (0)) 2 ), 
which implies ( t,p,p small) 

1 _ 2 1 _ 

(^ +ct )'V0| 2 ^ E X ab e a (0)e fe (0) < (_--rf)|V^| 2 , 



a,6=l 



t — U 



(2.33) 
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and (J23D yields 

2 



(XabU + X ab u)e a ((f))e b ((f)) 

a,6=l 

2 2 

= 2{t -u)J2 XaM^^) + ^ " m2 )^ H X afe e a (0)e 6 ((/>) (2.34) 

a, 6=1 a, 6=1 

2 

- 2m« k ab e a ((j))e b ((j)). 
0,6=1 

As A; a f, is bounded, combining (|2.33[) . (2.10) and lemma 2.4 we conclude 

2 

(2 + Ct)|V0| 2 < Yl (XabM + X ab ^a(^e b (ct>) < (2 - Ct)\ V</>| 2 . (2.35) 

a,6=l 

So we finish the proof of theorem 2.1. 



3 Non-concentration of the L 6 part of the energy 

In this section, we will prove lemma 1.2, which is essential to prove global existence and 
regularity. First we introduce some notations. 

Let zq = (0,0), be the vertices of the backward geodesic cone, then 

Q(z ) = {(t,x) £ [t ,0) x M 3 : u<0, t < 0}, 

denotes the backward geodesic cone, if to < s\ < S2 < 0, set 

Qll = Q(zo) n ([si,s 2 ]), 

and 

M tl = dQH = {(t,x) € Q s s \ : u = 0}, 
denotes the mantle associated with the truncated cone _ 

Q(s) ={i£l 3 : u < 0, t = s} 

denotes the spatial cross-sections of the backward cone Q(zq) when the time is s. 
Define the energy of problem (jl.ip 

Ex{t) = \l L 2 t +g ij (t,x)d i( j>d j( l)+^-) dx. (3.1) 


As we have showed in section 2 that d#u, VMt are bounded locally, then u t is continuous 
and together with ()2. 1|) we have 

lim m(t,x) = — = — = 1, (3-2) 

*,x-H) hm t , x -^o o t u(t,x) hm t , x ^o Otu{t,0) 
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that is: m = 1 + 0(t). So when t is small, Ei(t) has a equivalent form 

E(t) = -! (m' 1 (L(0)) 2 + m(L(0)) 2 + (m + m" 1 )^ 2 + - + J" - 6 ) eh;, (3.3) 
4 Jr3 v 3 / 

where |V0| 2 = (ei(0)) 2 + (e2(</>)) 2 , and (if = ^/\g\dx is the volume element corresponding 
to the metric g. Denoting the energy density 

e(t) = ^ (m- 1 ^^)) 2 + m(L(<A)) 2 + (m + m" 1 )^ 2 + "^"V ) . 
We then define the energy flux across M*: 

Flux 1 {(f>,M t s )= / -^-^ , ^ jy 3; = ~ ^-d^, (3.4) 

• / (^) 2 + E?=i(^) 2 



where g?z/ denotes the induced Lebesgue measure on M*. Similar to the energy, it has an 
equivalent form when t is small 

^M.)=f W±MMd, (3.5) 

where da = \f\g\dv denotes the volume element corresponding to the metric g on Mj , and 
it implies 

Fluxi&Ml) > 0. 

Lemma 3.1. When t is small, E\(i) and Fluxi(<p, M*) are equivalent to E(t) and Flux((j>, M*) re- 
spectively, that is: ±2 E(t), Flux i ((f), M*) ±2 Flux(<p, M*). 
Proof. Since 

L = m _1 (^ + 7V),L = m($ - TV), 



we get 



so 



And 



d t = -{m^L + mL), (3.6) 

W0 2 = [^K'iW + mLW)] 2 

= i(m^ 2 (L(0)) 2 ) + 2L(0)L(0) + m 2 (L(<A)) 2 ). 



<V</>, V0> = -{dttf + g^diWrt = ( ei (4>)) 2 + (e 2 (0)) 2 - L(0)L(0), 



which yield 

<fW^> = |V0| 2 - L(</,)L(4>) + (M) 2 , (3.7) 
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then we get 

Ei(t) = ' 
Fluxi(4>, Ml 



m- 2 (L(<j))y + m\L{<t>)Y + 2|V0|^ + -<f ) dx, 



Ml 



■((d t u(j) t ) 2 - 2d t u(t)tg %j diudj 



+ 



d t u 



iMW + lO- 28 t u 



' 'g ij dMf 



dv 



1 



** V(^a) 2 + Ej=i(^u) 2 

6 



f (^)) 2 + ^(IWI 2 - 2^)£(0) + (3^) 2 ) 



6m 2 V <9 t u 



di/ 



/ 

J Ml 



1 



m 
Y~2 



(^)) 2 + — (|V0| 
2m 



m-^a^ + N{4>))m{d t <l> - N{(j>)) + (^0) 2 ) + ^ - J- (A^/>)) 2 



3 



df, 



2V(^M) 2 + E-=i(^ 



together with (|3.2p . we obtain the result. 

To finish the proof, We shall require several other lemmas. The first is standard and 
says that the energy associated with our equation is bounded. 

Lemma 3.2. If G C°°([t ,0) x M 3 ) is a solution to (1.1), then Exit) or E(t) is bounded 
for all t < t < 0. Additionally, if t < s < t < 0, then 



Flux{(/>, Ml) -)• 0, when s, t ->■ 0. 



(3.8) 



Proof. To prove the boundedness of energy one multiplies both sides of the equation 
4>tt — g|r(<7 y (i 5 x )4>j) + <fr 5 = by dt<p to obtain the identity 



d_fcj> 2 +g^{t,x) 
dt 



+ 



1 



d t g l3 (t,x)(j)i(t)j 



_d_ 

dxi 



4> t g lj {t,x)4) 3 



0. 



(3.9) 



Thus, 

A2 i n ijU A$ 



Jr3 V J 6 / J R 3 2 1 ^ 



<j ><^(t,x)<^)dx = 0. 



9*. . 
Jrz dxi 

And since the last term is always zero, by the divergence theorem, due to the fact that 4>(t, x) 
for \x\ > C + t, (EHUD implies 

dtE^t) <CEx(t), 

which means 

Ex{t) < E x {t )e c ^\ 



15 



so Exit) or E(t) is bounded, as desired. 

To prove the other half of lemma 3.2, we integrate (|3.9p over and arrive at the "flux 
identity" : 



1 f tot, x) + g ij (t, x)di(f>(t, x)d j( f>(t, x) + X) ) dx + F/nxi (0, M l s ) 
Q(t) v a y 



2 



Q(s) 

d t g tj (T, x)di<f>(r, x)dj^(r, x)dxdr, 



__ 1 
~~ 2 

that is 

ft 



Ei((j), Q(t)) + Fluxi((j), M*) - Ei(</>, Q(s)) < C(t ) J EM, Q(r))dr, (3.11 
where C(to) is a constant depending on to- And it means 

Q(t)) - C(t ) /"* EM, Q(r))dr + Flm^, M s *) 
J t 

< EM, Q(s)) - C(t ) f EM, Q(r))dr, 

./tn 



(3.12) 



which implies Ei(<p, Q(t)) — C(to) // Ei(4>, Q(r))dr is a non-increasing function on [to, 0). It 
is also bounded as we have showed above, hence Ei((f>,Q(t)) — C(T) J t E%((j), Q(r))dr 
and E\{(f), Q(s)) — C(T) -Ei(</>, Q(r))dr in f)3.11j) must approach a common limit. This in 
turn gives the important fact that 

Fluxi(4>,M l s ) ->■ 0, when s,t ->■ 0, 

thanks to lemma 3.1, we complete the proof of lemma 3.2. 

To prove lemma 1.2, we need to introduce the energy- momentum tensor II as a sym- 
metric 2-tensor by 

U(X, Y) = X(4>)Y(4>) -^<X,Y> |V</>| 2 , 

n a/3 = d a (f)dp(j> - ^g a/3 \V(j)\ 2 , 

where X, Y are vector fields and <f> a fixed C 1 function. Then we have 

n(L,L) = {L(<p)f, II(L, L) = {L{4>)) 2 , 
H(L, e a ) = L(4>)e a (4>), U{L, e«) = L(4>)e a (4>), 

n(L, L) = L(0)L(0) - I < L, L > |V</>| 2 = + |V0| 2 = |V</»| 2 , (3 13) 

n(e a ,e 6 ) = e a (0)e 6 (0) - - < e a ,e fe > |V0| 2 

= e a (0)e 6 (£) - ^M|W>| 2 " L{<t>)L(<t>)), 

where d a b denotes the Kronecker delta function. 

We also need a key formula showed as a lemma below. 
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Lemma 3.3. Let be a C 1 function and II be the associated energy-momentum tensor. 
Let X be a vector field, and set P a = HapX", then 

divP = D a P a = a g <px(<f>) + hr^ x \ a p, (3.14) 
where O g is the wave operators associated to the given metric g and has formula as follows: 



^=\9\- l/2 d a { 9 a %\ 112 

1 1 (3-15) 

= -d tt <f> + di(gV(t, x)<t>j) + -g VJ g lm d m g l3 d^ - -g ij d t9ij d t cf>, 

where \g\ is the absolute value of the determinant of the matrix (<7 a «) and (g a/3 ) its inverse 
matrix. 

For the proof one can read pQ. 

We then construct a multiplier: ^(uL + uL) + 1, which is close to the Morawetz multi- 
plier tdt + rd r + 1, and setting Y = \{uL + uL). 

Following Christodoulou and Klainerman [3], the deformation tensor of a given vector 
field X is the symmetric 2-tensor ( x \r defined by 

Z) = tt(Y, Z) =< D Y X, Z > + < D Z X, Y > . 

In local coordinates 

^afi = DaXp + DpX a , 

as 

Vu = 2Vt - Vu = -2d t + L 

= -(m~ l L + mL) + L = -mr x L + (1 - m)L, 
then we can compute the deformation tensor of Y = \ (uL + uL) as follows 

(y) vrLL = 0, {Y \ LL = -2 + 2um, 

— — m 

^TLL = 4(1 - m) - 4uM, {Y) 7TLe a = (Va ~ %)u, 
{Y) KLe a = ZaU + 27/ a U, {Y) 7r ea e b = X al) U + XabU- 



Also 



divY = g afi < D a Y, ep >= g afi < D a ^(uL + uL), e p > 

1 . V _1 

= 2 ViaaU + Xbbu + x^u + X bb u) + 1 + m -uu. 



(3.16) 



Combining (|1.1|) and (|3.15|) . we get 



D 90 = <P 5 + ^ 9 lm dm9iM - \^Qt9iA^ ( 3 - 17 ) 
together with (|3.14p . substitueing X with Y we arrive at 

divP = D a p a = cvy($ + in^( y V Q/3 

= + \9 ij 9 lm d m 9iA<l> ~ l^dtg^Y^) + ~ir^ Y V Q/3 

j.6 -i -i I (3.18) 

= Y(y ) + (- 9 *V m a m «^ - -gfdtgijdtfiYte) + -ir^n^ 

= divtff) - ^divY + (^V™d m5ii ^ - \g ij d mj d^)Y^) + \TSf#V\ aP 
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where P a = H a gY^, and it means 

_ div{P _ l^Y) = ±6 6 divY - Y{6)^g lm d m g %3 d^ - \g K j d t9ij d t 6) 



\u^ Y \ a ^R(t,x). 



(3.19) 



By ([321) and dSJZD ; we have 
1 .. ,_.l 



° g (^ 2 ) = div{V(-0 2 )) = div(0V6) =< D a 6V6,d a > 



d a {4>) <vs,d a >+6u^ 

'0 -> -T^gf, 



g^d a {6) <V6,ds>+6n 



\V6\ 2 - + / + <t>^g ij g lm d m9ij di<t> - ^g ij d t9ij dt<t>), 



so 



-div^) = -\V6\ 2 + L(0)L(</>) - / - ct>(lg ij 9 lm dm9iA<t> ~ \g ij dt9iA<t>)- ( 3 - 2 °) 
Adding (I3J91) and ^3HD , we get 



diviP - -<p 6 Y + 0V0) = R(t, x) 
6 



- |V0| 2 + L(<£)L(0) - / - 0(^V m 3 m ^</, - ~g ij d mj d t <j>) =R(t,x). 



(3.21) 



Integrating the identity (|3.2ip over the truncated geodesic cone Q%, S <T < 0, we arrive 
at 

f l. f <P- \(h & Y + 6V 6, Vu > 
- <P- -6 6 Y + 6V6, -d t >dv- 6 y ^ v ^ 

Jq(t) $ Afj V(9 t u) 2 + WW 2 

+ I <P- -6 6 Y + 6V6, -d t >dv= [ R(t, x)dvdt, 
JQ(S) 6 J q t 



s 



that is 



/ n(y, < l<f> 6 Y - d t >dv- [ u(y, d t )- < I^y - 6Vcp, d t > dv 

JQ(T) J OAS) O 



'Q(T) u ^Q(S) 

J M$ V(dtu) 2 + (g ij diu) 2 Jqt 
By (|3.6p . we have 

U(Y,d t )- < ^6 6 Y-6V6,d t >=H(^(uL + uL), ^{mT 1 L + mL)) 

- < -6 6 -(uL + uL) - 6V6, -(m~ l L + mL) > 
6 2 2 

mii /r , ,,>2 w / T /,s\2 / "U mils.— ,.? . u mus , fi 

i 777 



(3.22) 



(3.23) 
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and 

< P - 1<P 6 Y + <f>V<f>, L >= U(Y, L)— < \<\?>Y - 0V0, L > 
6 6 

= IL(-(m~ 1 L + mL), L)- < -(() 6 -(uL + uL) - cj>V(/), L> (3.24) 

= \um 2 + \u{mf + ^ + <pl(4>), 

then ()3.22p becomes 



/ 



2m V W 2 ^- WJ 7 m t V(^) 2 + (5 ij '5iM) 2 



(3.25) 



- / ^(L( < A)) 2 + ^(L(^)) 2 +(^ + ^)|W| 2 + (^ + — W 6 
y Q(s) L 4 y - yV " 4m 1 Mm 4 Jl Vl V 12m 12 ;V 

+ — 4>L{^) + — <t>L(4>) ]dv = / R{t,x)dvdt, 
2m 2 Jq t s 

where Q(S) = {x E R 3 : u < 0, t = S}. Noting that u = on the mantle Mo, and 
when S, T is small enough we can let m = 1 for the error margin is nothing but 0(t 2 )E(t), then (13T2"5T) 
becomes a little simpler form 

/ [^(0)) 2 + f (M^)) 2 + |lW| 2 + + jW) + ^£fo)]di, 



MJ \/(^) 2 + (ff^u) 2 
JQ(S) 4 4 2 O 2 2 



(3.26) 



i?(t, z)dudi 



Denote 



i = J qt [jimy + jim) + ^<p\ 2 + ^ + + -mm^ 
u t(L(^)) 2 +0L(0) ^ = r t{mf + ^ 



mJ ^{dtu) 2 + (g lj diu) 2 Jm° s y/(dtu} 2 + (g^diu) 2 



m=- I [- A {mf + = {mf + |iwi 2 + 1 <t? + Um + \m<t>)]^, 



IQ(S) 

then (|3.26p becomes 



I + Ih- Ih + III = [R(t,x)dvdt. (3.27) 

JQ S 
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Let us estimate the right-hind side of (|3.27p first. 
rr«^(XL V „A3'tt (Y)„ 

= (wu-1- — )|W| 2 + (l-^-^)(^(0)) : 



m 

2 2 

- ~ %)uL{4>)e a {4>) + ^2 (Xabu + x a6 u)e o (^>)e 6 (0) 

a=l a,b=l 

1 _ _ (3-28) 

- ^(XaaM + + X aa ^ + X 6f) w)|V0| 

+ ^(XaaU + XbbU + x aa u + x bb n)L(0)L(0) 

2 

-^(CaM + 2 2o «)L(^)e a (^). 

a=l 

Combining (|3.16p (|3.19p (|3.2ip with (|3.28p . and set m = 1 (will not influence our result) we 
get 

R(t, X )dvdt = [ t ^^ + ^ + K aa - + X bb u) + 2-uu _ 2 } f_ 
Jqt 6 3 6 

+ (^(XaaU + + X^U + ~ 1 ~ ^ ~ 2 )) I V< ^l 2 

2 



2 

o,6=l 



+ (1 - j{Xa a u + XbbU + X^u + ^ h u))L{4>)^) - -wu(L(</))) 2 ^ 

1 2 2 

+ 2 " %)^ L ^) e M + J^fe + 2r/ o u)L(0)e a (0) 

a=l a=l 

--- dvdt. 



Also we have 



<C{T-S){ ^dvY{ du)3[(/ (d t 0) 2 d^ + (/ (6f^) 2 dv)s] 

JQ(S) ^Q(S) JQ(5) 

<C(T-S)|5|(£;(^Q(5)))5, (3 - 30) 

iT Z ^ _■ 

s 

<C\S\(T-S)(E(<f>,Q(S))). 
Combining (2.5), (2.6), (2.8), (2.9), (2.11), (2.12), (l23Tjl . (|3T29l) and (|330D . we get 

fl(t, x)dt;di < C|5| (T - S) Q(5))) + C(T - S)|S| Q(5))) 1 . (3.31) 
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On the surface Mj where u = 0, we have 

= t(m~ l dt4> — g^diudjcj)) 2 + (j)(m~ 1 dt4> ~ 9 K 'diudj4>) 

= — (u — t){<g % i diudjrf) — m~ l dt(p) 2 — <p{g %J 'diudj(j) — mT l dt(, 

If we parameterize M$ by 

y^(f(y),y), y€Q(S), 

then by u(f(y),y) = on Mg, we have 

Utfi + Ui = 0, 

fi = = -mui, 



and let ip(y) = <f>(f(y),y), then da = yj \d t u} 2 + (g^diu) 2 dy and 

= 4>tfj + 4>j, 

which implies 

g tj diudj^ = 4> t g lj diufj + g lj diudj(j) 

= —m(ptg l3 diudjU + g %3 diudj(f> 
= -m<p t (d t (f)) 2 + g tj diudj(f) 
= —mT 1 ^ + g lJ diudj<j). 

Thus, a calculation gives 

Ih = - f [(u - S) (g ij diudjip) 2 + ijg ij diudj^} dv 
JQ(S) 



r {(u-s)g i w m d j ^ + ^) 2 dv | r 

Jois) u — S Jc 



2 



+ tpg^diudjijjdv. 



IQ(S) U — O 7q( S ) u- S 



Integrating by parts we see 

/ tpg 13 diudjtpdv 
JQ(S) 



'Q(S) 

g^diudA^dv 
Q(S) 2 



[ [djil^dM 2 ) - l^djtfidiu)] dv 
JQ(S) 1 1 

[ fMm ^da-f U 2 d 3 ^ dl u)d V . 



Note that 



n g u = div(Vu) = -div(L) = ~x u - X 22 



-d a u + d j (g ij (t,x)d i u) + -g ij g lm d m g i:j diu- -g ij dtgijd t u, 



(3.32) 



(3.33) 
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which yields 

dj{g lJ (t,x)diu) = -x u -X 22 + d a u- g lm d m gijdiu + ^g l:l d t gijd t u, 
then from (2.6) and (2.7), we have 

2 +Ct + C<d, (g ij diu) < — Ct + C. (3.34) 



u — t u — t 

Combining (13321 . (13331 and (1331]) we get 

Ih = - [ (fe-W^ + tf dv + f ^Mm „,, 2J „ 

-(CS + C) [ ^ 2 dv 
Jq(s) 

(u-S){- m- 1 ^ + gVdiudrf + ^) 2 + (CS + C)^ 



da 

MO 

2^£-=i(^) 2 

= /■ g(L(0) + |) 2 + (cg + c)0 2 dj | r g^g^ 

A/g V(^) 2 + (g^diu) 2 JdQ(s) 2^Y, 3 j=i( d m) 2 

<C\S\ [ {mfda + C [ (t4| +1+ I^D^d^ 



<C\S\Flux(cj),M s )+C(\S\ + \S\ 2 + \S\ 3 )( ( do) 



Ml 



2^E J 3 =1 (5 J n) 2 



For JiT, a computation gives (also let m = 1 



2 A /E =i(<%) 2 



^W>)) 2 + =(m) 2 + f |V0| 2 + f / + + ^L(0) 

~ ' " ~ l d t 4> - g lj diudj(p) 2 + j (md t 4> + m 2 g lj diudj(f) 2 



(3.35) 



4 ^ m 

+ f [ - (d t <P) 2 + g ij d^djcj) + (m^dtcj) - g^diud^) (md t cP + m?g ij diudrf)] 
b 

= § + <f^#^ + y ) + & + (« - S)^^^) . 
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For the second term on the right-hand side, using Cauchy- Schwartz inequality we have 
<t> t {(t>+(u-S)g ij d i ud j (t>) 



-'^T + 2|5p 

+ (u- S)gV diudjcf) 



2 



2 



<\s\[^ + 



2 ' 2(u-S) 2 J 

^, ol ^ , |5| r <£ 2 , 4i _ . iia Ma x , l^djud^ 

~ ' '~2~ ~2~ >-ff) 2 5 J 5f W + — u _ s \ 

^■ g i^ , |g| r ^ , -2 , , ^djud^ 



As 5 < 0, we get 



| (0? + cf^a^ + ^) + 4>t {<t> + (n - S)g^diud j( j>) 
~ 6 ~ 2(u-S) 2 u-S ' 



so 

A6 



/// = -/ [| {4 + <^0fi^ + ^-) + <k (<f> + (u - S)g ii d i ud j( f)}dv 
JQ(S) 1 6 



'Q(S) 6 V2 JQ(S) [U- b) z Jq {S ) u-s 

Together with ()3.34p . a similar computation gives 



(f)g tj diudj (t> 
Q(5) M - 5 1 



dt> 



g lJ diudj(^ 



-dv 



x u-S J . o<ff\ 2 V« - S. 



***** -du-l ^-A^- + CS + C)d v 



80(5) -2S a /eU(%I) 2 2 (^ " 5 ) M " S 



+ J Q(S) 2(u-sy dv - 

Combining (|3.36|) and (|3.37|) . we get 



(3.36) 



Q(5) v Ml- J ' JQ(S) 

f (f 1 ^ d ud ucjP' f (jp 1 ■ ■ 

f (t> 2 g lj djudju ^ 
0(5) 2(u-Sf 



III > \S\ [ £dv- f tiMML da - (CS + CS 2 ) / -^d„. (3 38) 
~ J Q{ S) 6 J dQ{s) 2^~Jp~uf JQ(S) 2(M - S) ( 3 - 38 ) 
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Using Holder's inequality it is easy to see that 
I -- 



IQ(T) 

<C\T\E(<f>,Q(T)) + C\T\ 



2 

fdv 



Q(T) 



Q(T) 



< C\T\E(^Q(T)) + C\T\E{^ Q(T))s 

>2 



Ih 



t(L(d>)) +6LU) , , , nx ■ , , n, 

/ o d i , =k d<T ^ C\T\Flux{^ M°) + C\T\Flux{<j>, M°) 



(3.39) 



Now, we combine (13371) . (133Tj) . (13351) . (13351) and (^391) to obtain 



|5| 



b -dv<III + f 3 ij djudM 2 



Q(S) 



8Q(S) 2E,=i(3j«) 



da + (OS + OS 2 



R(t,x)dvdt + 



-du 



= -/ - Ih + 7/ 2 + 

+ (CS + CS 2 )[ 

Jq(5) 2(n - 6) 

<c|r|(£;(0,Q(r) + ^,Q(r))§) 

+ C\S\ (Flux((f), Mg) + Flux(cj), Mg)^) 
+ C\T\ (Flux(cj), My) + M^)f ) 

+ C|S| (T - 5) (£(0, Q(5)) + (E(0, Q(S))) § ) 
+ (C,S 2 + C,S 3 ) (£(</>, Q(S)))^, 



Q(S) 2(u- S 1 ) 
da 



dv 



g lJ diudjucf) 2 



9Q(S)2E-=ife) 2 



and then the result of lemma 1.2 follows as we can choose T 



-S 2 . 
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